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A series  of  l-D  calculations  simulating  the  impact  of  the  flyer  plate  on  the 
granite  block  in  the  HUSKY  PUP  experiment  was  made  to  evaluate  the  effects  of 
c(|uation  of  state  (EOS)  modeling,  energy  levels  and  distributions,  and  numeri- 
cal parameters  on  the  resulting  time  of  arrival  (TOA)  of  the  shock  wave  at 
stations  in  the  block.  This  objective  was  to  provide  a basis  for  judging  the 
significance  of  comparisons  between  TOA  vs  depth  measurements  and  results  of 
more  complex  2-1)  pretest  calculations  of  the  event.  The  total  energy  in  the  _ 
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INTRODUCTION 


1 . 1 Background 

Tlie  HUSKY  PUP  underground  nuclear  test  included  an  experi- 
ment in  which  a curved  flyer  plate  consisting  of  a sandwich  of 
high  and  low  density  materials  was  accelerated  by  the  nuclear 
explosion  and  impacted  onto  the  face  of  an  8- ft  diam-'ter  cylin- 
der of  granite  located  30  cm  from  the  center  of  the  source 
(figure  1 ) . 

Prior  to  the  test,  2-D  rad i at  ion  - hydro  code  calculations 
of  the  event  were  performed*,  with  the  intent  of  evaluating 
the  ability  of  current  numerical  techniques  to  predict  the 
dynamics  of  acceleration  of  bomb  fragments  or  debris,  and 
their  subsequent  impact  on  the  ground  in  a near-surface  nuclear 
burst.  The  basis  of  tliis  evaluation  was  to  be  comparisons  of 
calculated  characteristics  of  tlie  sliock  wave  in  the  granite 
block  (times  of  arrival,  peak  pressure,  and  waveforms  vs  dis- 
tance into  the  block)  with  measurements  of  these  (juantitics 
obtained  from  the  HUSKY  I’UI’  experiment. 

I'igure  2 shows  the  comparison  of  the  calculated  and  mea- 
sured times  of  arrival  (TOA)  vs  depth  in  the  granite  block. 
Unfortunately,  the  attempts  to  measure  peak  pressures  and 
waveforms  wore  not  successful.  However,  the  attenuation  of 
peak  pressure  with  depth  in  the  block  can  be  inferred  from  the 
experimental  TOA  data  and  from  the  equation  of  state  models. 

This  interpretation  of  the  data  is  given  in  Appendix  A. 

1.2  Objectives  of  the  Current  Study 

Despite  some  fairly  substantial  differences  between  the 
2-D  hydrocode  tecluii((uos  and  models  (Lagrangian  vs  Eulcrian  grids, 

* These  calculations  were  made  liy  Lilley  at  LASlI  , BJi  i Ip)' 

Systems,  Science,  and  Software^,  and  Schlaug  at  Science 
Applications,  Inc.-^. 
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igure  1.  Schematic  .of  HUSKY  PUP  Granite  Block  Experiment  and  Geometry  of 
2-D  Calculationsl>2>3 


T&  (Hia^fecci  "*  “** 


Figure  2.  Experimental  Time  of  Arrival  Data  Compared  with  the 
Values  Obtained  In  Preshot  2-D  Calcvilatlons . 
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I’yatt  vs  LASL  EOS,  for  example),  two  of  the  calculations  (S'^ 
and  LASl.)  gave  TOA  vs  deptli  curves  which  fall  reasonably  close 
to  the  HUSKY  PUP  measurements.  This  could  he  for  two  reasons: 

1.  The  2-D  calculated  results  were  fortuitous 
(.  i . c . , one  or  more  of  tlie  calculators  made 
compensating  errors). 

2.  The  TOA  vs  deptli  in  the  granite  block  is 
insensitive  to  major  differences  in  tlie  basic 
assumptions  and  ca  1 cu  1 a t i ona  1 tecliniques. 

Tliese  poss  ill  1 i t i es  , plus  tlie  fact  tliat  tliere  are  some  con- 
sistent discrepancies  in  I'igure  2 lietween  tlie  results  of  the 
various  ca  1 eu  1 a t i ons  , as  well  as  between  tlie  calculations  and  tlie  mea- 
surements, indicate  a need  for  furtlier  evaluation  before  any 
conclusions  are  drawn  about  the  validit)’  of  tlie  2-D  code  ca  1 - 
cul a t ions. 

One  way  to  make  tliis  evaluation  would  be  to  repeat  the  2-D 
calculations  several  times  so  as  to  assess  tlieir  sensitivity 
to  various  input  parameters.  However,  because  2-11  calculations 
are  exjiensive  and  complex,  it  was  decided  to  use  a scries  of 
1-1)  calculations  for  tlie  sensitivity  stud>’  which  is  reported 
herein.  The  specific  objective  of  this  study  has  been  to 
evaluate  the  sens  i t i \’ i t >•  of  calculated  times  of  arrival  to 
uncertainties  or  variations  in  a number  of  parameters  involving 

the  pro|)erties  of  the  flyer  jilate  when  it  impacts  the  granite,  the 
high  pressure  propei'ties  ascribed  to  the  granite,  and  the  numerical 
parameters.  Irom  this  evaluation  of  TOA  sens  i t i v i t\’ , it  becomes 
possible  to  judge  the  significance  (or  lack  thereof)  of 
agreement  between  calculated  and  measured  times  of  arrival  in 
t he  granite  block. 


t h 
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1 . 3 Approach 


The  calculations  were  all  made  using  CRALE-1,  a 1-D 
arbitrary-Lagrangian-Euler ian  finite  difference  code.  The 
test  geometry  was  modeled  as  a 1-D  spherical  system  with  the 
granite  block  face  at  a radius  of  50  cm. 


The  flyer  plate  was  assumed  to  be  a uniform  material.  All 
of  the  calculations  were  started  at  the  instant  of  flyer  plate 
impact  on  the  granite;  flyer  plate  acceleration  by  the  nuclear 
burst  was  not  treated. 

.\  baseline  set  of  conditions  was  first  analyzed  which 
gave  a TOA  vs  depth  curve  that  follows  tlie  e.xper imenta  1 data  with 
good ' f i dc  1 i t y , es[ieciall>’  considering  the  1-1)  nature  of  the  analysis 
(.See  Section  2).  for  each  of  tlie  other  cases  analyzed,  the  TOA 
vs  depth  curve  was  compared  against  the  baseline  curve.  The 
difference,  expressed  as  a sensitivity  factor,  was  used  as  a 
measure  of  tlie  change  in  the  calculated  results  due  to  the  para- 
meter under  consideration.  Talile  1 summarizes  tlie  principal 
variables  which  were  incorporated  into  the  study.  A synopsis 
(originally  prepared  by  the  first  author  while  at  RDAl  of  the 
eight  EOS  models  investigated  in  this  study  appears  in  Appendix  B. 
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TABLE  1 

SUMMARY  OF  PARAMETERS  IN  THE  1-D  CALCULATIONS 

• FLYER  PLATE  VARIABLES 
Uniform  Imjiact  Velocity 

(40,  45,  47. 5,  50,  60  cm/sec) 

Velocity  Gradients  in  Plate 

(40  to  60,  55  to  65,  -40  to  70) 

Plate  Deiisity 

(1.325,  2.65,  5.3  g/cm^) 

Internal  Energy 

(~0,  100,  300,  500,  1450  Mh-cc/gm) 
Equation  of  State  for  Plate 
(Pyatt,  Schuster,  Tillotson) 

• EQUATION  OF  STATE  MODELS  FOR  GRANITE 

Pyatt 

Green- 1 , - 2 , - 3 
Kal itkin 
LAST 

Schuster 
Ti 1 lot  son 

• OTHER  GRANITE  BLOCK  AND  EOS  VARIABLIiS 

Calculated  Pressure  Level 
Release  Paths 
Preheating  by  Radiation 
EOS  Interpolation  Method 

• NUMERICAL  VARIABLES 
Zon i ng 

Artificial  Viscosity 
Energy  Iteration  Method 
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1.4  Conclusions 


The  major  conclusions  are: 

1.  The  most  important  parameter  in  the  calculations 
of  the  HUSKY  PUP  granite  block  experiment  (and 
presumably,  therefore,  in  debris  slap  from  a 
nuclear  burst)  is  the  total  energy  of  the  im- 
pacting mass.  The  distribution  between  kinetic 
and  internal  is  of  secondary  importance;  internal 
energy  is  about  3/4  as  effective  as  kinetic  in 
driving  the  initial  shock  into  the  block. 

2.  Over  a wide  variation  of  material  models,  the 
maximum  differences  in  the  close-in  (P  > 1 Mb) 
arrival  times  due  to  the  various  EOS  models  for 
the  flyer  plate  and  granite  are  equivalent  to  a 
change  of  only  301  in  the  initial  energy  of  the 
plate.  Since  the  resultant  motions  in  the 
granite  block  should  scale  approximately  as  the 
cube-root  of  the  energy,  this  difference  is  small. 
Hence,  the  simplest  reasonable  material  model  should 
suffice  to  calculate  the  high  pressure  properties. 

3.  The  difficulties  in  obtaining  experimental  data 
(especially  waveforms)  in  the  extreme  pressure 
regime  of  the  granite  block  in  the  HUSKY  PUP  test, 
and  the  apparent  insensitivity  of  calculational 
results  to  input  assumptions  in  this  regime, 
suggest  that  future  efforts  for  code  validation 
should  emphasize  comparisons  with  test  data  taken 
in  lower  pressure  (below  1 Mb)  regimes,  where 
material  properties  may  be  more  significant. 

r 
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If  the  correlation  between  the  TOA  sensitivity  factor,  | 

I,  (as  defined  in  Section  3)  and  total  energy  obtained  in  this  | 

1 -n  study  is  extended  to  the  2-0  calculations,  an  additional  | 

conclusion  may  be  inferred:  | 

4.  The  LASL,  , and  SAI  calculations  differ  | 

from  the  experiment  by  18,  23,  and  64%  in  j 

effective  energy.  Thus  the  LASL  and  S^  results  | 

are  probably  within  the  uncertainty  limits  of  i 

the  HUSKY  PUP  source  yield.  Further  refinements  5 

of  those  calculations,  in  order  to  better  match  j 

the  experimental  data,  do  not,  therefore  appear  j 

to  be  necessary  or  fruitful.  We  conclude  that  I 

the  technology  represented  by  these  2-0  calcu-  j 

lations  is  adequate  to  predict  debris  slap  I 

motions  in  cratering  and  ground  motion  analyses.  j 

The  substantially  lower  equivalent  energy  in  | 

the  SAI  2-0  results  is  not  readily  explainable  * 

by  differences  in  models,  techniques,  or  assump- 
tions. This  suggests  the  possibility  of  a 
calculational  error  of  some  sort. 


SECTION  2 


EASEL  INI-;  CALCULATION 

The  baseline  set  of  conditions  for  tiie  1-D  analyses  was 
as  follows: 

EOS  (for  both  flyer  plate  and  granite  block)  - Pyatt 

Elyer  plate  material  - granite 

Initial  density  (for  both  flyer  plate  and  granite 
block)  = 2.b5  g/cm^ 

1-lyer  plate  velocity  = 50  cm/nsec  (uniform) 

Elyer  plate  thickness  = 2.0  cm 

Internal  energy  - flyer  plate  = .005  Mb-cc/gm* 

granite  block  = .0001  Mb-cc/gm 

The  thickness  of  the  flyer  was  chosen  such  that  its  areal 
density  would  be  tiie  same  as  in  the  concentric  flyer  in  the 
experiment.  The  uniform  initial  velocity  of  tlie  flyer  was 
chosen  to  give  a total  kinetic  energy  rouglily  equivalent  to 
that  seen  in  ])rcliminary  calculations  of  tlie  experiment  by  S^. 

The  fl\-er  was  resolved  by  5 cells  of  equal  thickness 
(0.4  cm).  Lhe  cells  in  the  granite  block  increased  by  0.25°d 
per  cell  from  an  initial  thickness  of  0.4  cm  (at  the  flyer- 
block  interface).  At  a radius  of  210  cm  (180  cm  into  the  block, 
tlie  end  of  tlie  region  of  interest),  the  cells  were  about  0.9  cm 
thick. 

The  iOA  vs  depth  curve  obtained  using  the  baseline  condi- 
tions is  shown  in  1 igure  5**,  along  with  the  experimental  data 
and  with  the  results  of  the  2-U  calculations.  While  the  baseline 


* 1 Mb-cc/gm  = 10‘‘'  ergs/gm 

**  Incremental  methods,  like  finite  difference  codes,  do  not 
provide  a smooth  TOA  curve,  but  ratber  a series  of  discrete 
points  corresponding  to  the  shock  arrival  at  each  cell.  flic 
smooth  curves  shown  throughout  this  report  are  least-square 
fits  of  the  form  P = at*’  of  the  discrete  points  calculated 
between  bO  cm  and  180  cm  depth  in  the  granite. 


1 1 


3.  Comparison  of  the  Times  of  Arrival  for  the  1-D  Baseline 
Case,  the  2-D  Pre-shot  Calculations  and  the  HUSKY  PUP 
Experiment . 
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calculation  was  not  intended  to  exactly  match  the  test  condi- 
tions, it  was  gratifying  to  obtain  the  good  correspondence 
shown  in  Figure  3.  This  correspondence  helps  to  establish 
the  relevance  of  the  baseline  conditions  to  the  objectives  of 
the  study  and  also  sliows  that  1-D  simulation  of  the  experiment  is 
quite  realistic.  There  is  an  increasing  discrepancy  between 
the  calculated  and  measured  TOA,  with  the  measured  TOA  being 
about  10b  later  than  the  calculated  TOA  by  180  cm  depth.  This 
presumably  reflects  the  2-0  attenuation  effects  (e.g.  lateral 
relief)  which  would  be  expected  to  appear  at  late  times  in  the 
experiment . 

Figure  4 shows  the  stress,  velocity,  and  displacement  time 
histories  calculated  in  the  Baseline  Case  at  a representative 
depth  (95  cm)  in  the  granite  block.  Figure  5 is  a plot  of  the 
peak  pressure  vs  depth.  The  knee  at  about  60  cm  depth  is  caused 
by  a strong  second  pulse.  This  pulse  arises  when  the  rear  of 
the  flyer  plate  converges  and  rebounds  from  the  center  of  the 
1-0  splierically  symmetric  geometry  of  the  problem.  Such  a 
secondary  pulse  may  not  be  present  in  the  actual  test  geometry. 


D isp Laiemvitit  (cm)  Velocity  (cra/psec)  Pressure  (Mb) 

60.  too.  -.00  .60  1.00  0*00  2.50  6.00 
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Figure  4.  Stress,  Velocity,  and  Displacement  Histories  in  Baseline  Case 
at  95  cm  Depth  in  Granite  Block. 
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SECTION  3 


BASES  FOR  JUDGING  EFFECTS  OF  PARAMETRIC  CHANGES 

Where  single-parameter  variations  are  made  in  a sensitivity 
stud'  . and  where  the  meaningful  output  consists  of  a single 
quantity,  a common  procedure  is  to  determine  the  fractional 
change  in  the  output  quantity  as  a f unct ion  of  fractional 
changes  in  the  input  parameters.  This  approach  was  not  practi- 
cal in  the  present  study,  inasmuch  as  the  TOA  vs  depth  curve 
is  not  a simple  output  quantity,  and  because  some  of  the  input 
variables  involved  combinations  of  parameters  and/or  overall 
model  changes.  For  this  reason,  the  criteria  described  below 
were  used  for  judging  the  effects  of  input  variables  on  the 
results  of  the  1-D  calculations. 

3. 1 Time-of-Arrival  vs  Depth  Sensitivity  Factor  (Z) 

The  TOA  vs  depth  curves  for  each  case  considered  were 
compared  against  the  curve  for  the  baseline  case,  and  a non- 
dimensional  TOA  sensitivity  factor,  E,  was  calculated.  E is 
the  integral,  illustrated  in  Figure  , of  the  pevacntage  dif- 
ference in  arrival  time  between  the  two  cases  in  the  interval 
between  60  and  180  cm  depth,  normalized  by  that  interval  (i.e., 
120  cm) . 

E is  considered  as  positive  for  cases  where  tlie  TOA  is 
always  earlier  than  in  the  Baseline  Case,  and  negative  if  always 
later.  l-'or  the  few  cases  where  the  TOA  vs  depth  curve  crossed 
the  baseline  curve  between  60  and  180  cm  depth,  E was  computed 
as  the  sum  cf  the  absolute  values  of  the  partial  integrals  on 
each  side  of  the  crossing  point,  witli  tlie  sign  of  the  larger. 


in 


3.2  liffcctivc  Lquivalent  Kinetic  Energy,  E^ff,  and 

Effective  Energy  Ratio, 

Since  the  shock  propagation  velocity  is  directly  deter- 
mined by  the  impact  velocity,  it  is  obvious  that  the  impact 
velocity  of  the  flyer  plate,  and  hence  its  kinetic  energy,  has 
a major  effect  on  TOA  vs  depth  in  the  granite  block.  For  this 
reason,  it  was  decided  to  relate  results  of  those  cases  in- 
volving other  variables  to  the  results  obtained  when  only  the 
uniform  impact  velocity  of  the  flyer  plate  was  changed.  In 
this  way,  an  effective  equivalent  kinetic  energy,  Egfp.  could 
be  obtained,  as  well  as  the  ratio  of  this  value  to  the  total 
actual  energy  in  the  flyer  plate,  Egj-f/E^g^. 

To  do  this,  tlie  relationship  between  the  TOA  sensitivity 
factor,  Z,  and  the  energ)'  (essentially  all  kinetic)  of  uniform- 
velocity  flyer  plates  was  first  derived  by  fitting  a quadratic 
to  the  values  of  I and  for  the  Baseline  Case  (V'=50  cm/psec). 

Case  PI  (V=40)  , and  Case  l’2  (V'=60)  , i.e.. 


Z =-.3227  . 1.119  E^^^  - .7964 


(1) 


Eejn.  1 will  he  referred  to  as  the  standard  E-E^^^  curve  for 
un i form- ve 1 oc i t y , al 1 -k inet ic -energy  flyer  plates. 

The  value  of  Eg(-{-  for  each  subsequent  run  was  obtained  by 
inserting  the  value  of  tlie  TOA  sensitivity  factor,  Z,  calculated 
for  each  run  into  F.c|n.  1 and  solving  for  E^j-^  is  thus  the 

kinetic  energy  of  a fl)'er  plate  of  uniform  velocity  which  will 
produce  tlie  same  TOA  sensitivity  as  the  variable  or  variables 
under  consideration  in  a specific  l-D  case. 
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Figure  7 graphically  illustrates  the  determination  and  use 
of  Ejjff  This  is  a plot  of  the  TOA  sensitivity  factor,  Z,  vs 
total  energy,  (normalized  by  the  total  energy  of  the  Base- 

line Case).  The  solid  curve  in  Eqn . 1,  represents  the  standard 
relationship  between  Z and  E^^^  when  the  plate  velocity  is  uni- 
form and  virtually  all  energy  in  the  plate  is  initially  kinetic. 
As  an  e.xample  of  how  results  of  1-D  cases  arc  related  to  Eqn.  1, 
the  calculated  value  of  Z,  and  E^^^  arc  plotted  for  Case  P7. 

Its  total  energy,  is  79.1-6  of  the  Baseline  Case,  but  it 

produces  the  same  effect  on  the  TOA  sensitivity  factor,  Z,  as  a 
un  i form- ve  loc  i t y plate  having  72.3"6  of  the  energy  of  the  Base- 
line Case.  Thus  the  effective  equivalent  kinetic  cnerg)',  E^^^ 
for  the  F’7  flyer  plate  is  .725,  and  its  effective  energy  ratio, 
i ^cff^^tot’  .723/. 791  = .915.  E^^^/E^,^^  is  seen  graphically 

f in  Figure  7 to  he  merely  a measure  of  the  horizontal  distance 

I between  the  plotted  value  of  Z and  Ej.^^  for  any  case  and  the 

standard  curve  (lain.  1). 

I 
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SECTION  4 


RESULTS  OF  1-D  CALCULATIONS 

The  38  one-dimensional  (spherical  symmetry)  cases  analyzed 
in  this  study  are  listed  in  Table  2.  They  are  grouped  into  four 
general  categories  according  to  whether  the  parameter  being 
examined  was  associated  primarily  with  the  flyer  plate  (the 
P-cases)  , granite  equation  of  state  model  (E)  , other  granite 
block  and  EOS  variables  (B)  , or  computational  technique  (C) . 

Table  2 summarizes  the  most  important  results  of  the  study, 
including  the  TOA  sensitivity  factor  (Z)  , the  equivalent  effec- 
tive kinetic  energy,  E^^£,  and  the  effective  energy  ratio, 
E^ff/Etot>  defined  in  Section  3. 

4 . 1 Flyer  Plate  Parameters 
4.1.1  Impact  Velocity 

Initially,  the  simplest  and  most  obvious  parameter  to  vary 
was  the  impact  velocity  of  the  flyer  plate  on  the  granite  block. 

In  cases  BL , PI,  P2 , P3,  and  P4 , the  velocities  within  the  plates 
were  uniform.  These  velocities  were  varied  from  40  cm/qsec  to 
60  cm/psec.  The  TOA  vs  depth  for  these  cases  is  shown  in  Figure  8. 
As  expected,  increasing  the  impact  velocity  produces  a consistent 
decrease  in  the  TOA  at  a given  depth. 

Since  the  actual  flyer  plate  velocity  was  probably  not 
uniform  through  the  plate,  three  problems  (Cases  P5,  P6 , and  P7) 
were  run  in  which  the  velocity  increased  linearly  with  radius 
within  the  plate.  Thus  in  Case  P5,  the  front  (impacting)  surface 
of  the  plate  was  traveling  at  6.S  cm/pscc,  while  the  rear  was  at 
35  cm/psec,  with  linear  variation  in  between.  In  Case  P7,  the 
front  surface  was  at  70cm/ps'c,  while  the  rear  surface  was 
moving  at  -40  cm/pscc  (i.e.,  it  was  converging  back  towards  the 
source).  This  roughly  approximates  a velocity  profile  calculated 


TABLE  2.  SUMMARY  OF  1-D  CASES  AND  RESULTS 


by  S''  in  a preliminary  2-D  study.  TOA  results  of  these  cases 
are  compared  with  the  Baseline  Case  in  Figure  9. 

Figure  1 0 shows  the  results  of  all  of  these  cases  on  the 
Z-Etot  plane,  and  compares  them  with  the  standard  relation- 

ship (Eqn.  1).  The  three  non-uniform  velocity  cases,  P5,  P6, 
and  P7,  fall  only  slightly  below  the  curve  (i.e.,  the  TOA  arc 
generally  later)  indicating  that  when  the  kinetic  energy  is 
non-uni formly  distributed  in  the  flyer  plate,  its  effect  on  the 
TOA  in  the  granite  plate  is  less  by  a small  amount  than  an 
equivalent  energy  uniformly  distributed  in  the  plate.  This  is 
somewhat  surprising  in  Case  P7,  since  part  of  the  kinetic  energy 
in  that  case  was  directed  away  from  the  granite  block.  Even  so, 
the  effective  energy  ratio,  Eg££-/E^^^,  for  Case  P7  is  .915,  which 
says  tluit  the  strongly  non-uniform  nature  of  the  P7  flyer  plate 
resulted  in  a loss  of  effectiveness  of  only  8.5?>,  as  compared  to 
a uni  form- ve loc i ty  flyer  plate.  This  suggests  that  the  total 
kinetic  energy  in  the  flyer  plate  has  a much  more  important 
effect  on  the  TOA  vs  depth  in  the  granite  block  than  the  detailed 
distribution  of  the  kinetic  energy  within  the  flyer  plate, 

4.1,2  Density  and  Flyer  Plate  Equation  of  State 

The  second  flyer  plate  parameter  studied  was  its  density. 

The  pre-test  plate  consisted  of  sandwich  layers  of  high  and  low 
density  materials.  At  the  time  of  impact  with  the  granite  block 
the  actual  flyer  plate  had  presumably  undergone  several  stages  of 
compression  and  rarefaction  due  to  the  force  of  the  explosion, 
and  its  density  profile  was  clearly  uncertain.  In  the  Baseline 
Case,  the  flyer  plate  was  granite  with  a uniform  density  of 

= 2.65  g/cm-^.  In  Cases  P8  and  P9 , was  increased  and  decreased 
by  a factor  of  2,  to  5.3  and  1.325  g/cm^.  The  effects  of  these 
changes  on  TOA  vs  depth  in  the  granite  block  arc  seen  in  Figure  11 
to  be  quite  sma 1 1 . 
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Figure  11.  Variations  in  Time  of  Arrival  Due  to  Changes 

in  the  Initial  Density,  p^,  of  the  Flyer  i’late. 
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Similarly,  when  the  EOS  model  for  the  flyer  plate  was 
changed  from  the  Pyatt  model  (used  in  the  Baseline  Case)  to 
the  Tillotson  model  (Case  P19)  , the  change  in  TOA  vs  depth  was 
negligible.  Even  when  the  flyer  plate  was  changed  to  expanded 
Tungsten  (p^  = 2.65),  and  its  internal  energy,  IEq,  was  raised 
to  100  Mb-cc/gm,  (Case  P20)  , the  effects  on  TOA  vs  depth  in  the 
granite  block  were  small.  These  flyer  plate  EOS  comparisons 
are  shown  in  Figure  12. 

When  the  TOA  sensitivity  factor,  E,  and  the  normalized 
total  energy,  for  the  flyer  plate  density  and  EOS  cases 

are  plotted  (Figure  1.5),  it  is  seen  that  they  deviate  only 
slightly  from  the  standard  E vs  E^^^  relationship  (Eqn.  1).  The 
largest  deviation  is  for  the  low  density  flyer  plate  (Case  P9 , 
with  p=1.325).  For  this  case  the  effective  energy  ratio  is 
still  .978,  or  just  2.2°6  less  effective  than  a flyer  plate  of 
density,  = 2.65  gm/cm^. 

4.1.3  Internal  Energy  in  Flyer  Plate 

Tavo  sets  of  cases  were  run  to  investigate  the  effects  of 
the  internal  energy  of  the  flyer  plate  on  the  TOA  in  the  granite 
block.  In  the  Baseline  Case,  the  initial  internal  energy,  lE^, 
was  essentially  zero  (.005  Mb-cc/gm).  In  Cases  Pll  and  P12, 
lE^  was  increased  to  100  and  500  Mb-cc/gm,  respectively.  Since 
the  velocity  was  held  at  Vq  = 50  cm/iisec,  these  changes  in  inter- 
nal energy  raised  the  total  energy  in  Cases  Pll  and  P12  by  8.41 
and  42.2"6.  As  might  be  expected,  increasing  the  flyer  plate 
internal  energy  shifted  the  TOA  toward  earlier  arrivals,  as  seen 
in  Figure  14.  However,  when  the  TOA  sensitivity  factor,  E,  is 
plotted  vs  normalized  total  energy,  E^^^  in  Figure  15,  it  is 
seen  that  Cases  Pll  and  P12  fall  below  the  standard  E vs  E^^^ 
curve  (Eqn.  1)  for  a 1 1 - k i net i c - energy  cases.  This  indicates 
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TOA  Sensitivity  Factor 


th;it  when  internal  enci'gy  is  adJod  to  the  flyer  plate,  the  shock  ^ 

strength  produced  in  t!ie  granite  is  lower,  (and  hence  the  I'OA 

is  later)  than  would  he  the  case  if  an  equivalent  amount  of 

kinetia  energy  were  to  be  added.  More  specifically,  in  Cases 

I’ll  and  P12,  the  normalized  total  energy  was  raised,  through 

addition  of  internal  energy,  to  1.084  and  1.422,  respectively. 

The  effective  eiiuivalent  kinetic  energies  for  tliese  cases  was 
1.051  and  1.5(Hi.  Thus  increasing  internal  eyiergv  by  8.4-o  and 
42. 2o  produced  the  same  effect  as  would  be  achieved  by  increasing 
kinetic  energy  by  5.1o  and  50. (i?..  I’ut  another  way,  internal 
energy  in  these  cases  was  onlv  about  5/4  as  effective  as  an 
equivalent  amount  of  kinetic  energy. 

fo  sec  it'  this  effect  is  independent  of  the  assumed  initial 
kinetic  energy  level,  three  additional  cases  were  I'un  in  which 
tlie  same  -4(i  to  ”0  cm/usec  ramji  velocity  distribution  as  used 
in  Case  1'^  was  specified.  Internal  energy  densities  varied 
from  essential  tero  (Casr  T~)  to  100  Mb-cc/gm  (Case  lA),  500 
.'Ib-cc/gm  (I’lS),  and  500  .'Ib-cc/gm  (l’16).  fhe  fO.A  vs  depth  curves 
for  these  four  eases  are  eompaied  in  figure  1 (i  the  lOA  sensi- 
tixitv  factor,  Z,  and  total  energies  are  included  in  fii'.ure  15 
ihese  cases  are  nomin.illv  fit  b\'  the  dashed  line  in  figure  15,  wliich 
has  a somewhat  sm.iller  slojie  than  the  standard  curve.  This  confirms 
the  earlier  finding  that  addition  of  internal  energy  to  the  fl\’ei' 
pl.ite  has  .1  smaller  el  feet  on  tlian  tlie  addition  of  i in-  s.ino 

amount  of  kinetic  energy.  for  leases  I'lA,  IMS,  and  1’ 1 b , the  I 

ch.inge  in  I'O.V  sensilixit)-  factor,  T,  was  again  equivalent  to 

that  e.xpected  if  the  added  internal  energv  had  been  replaced  by  I 

about  5/1  as  much  kinetic  energy.  | 

Die  indication  that  internal  energv  in  the  fiver  plate  is  j 

equivalent  to  about  5/4  as  much  kinetic  energv  was  seen  even  | 

when  the  initial  velocitv  was  set  to  zero,  and  the  internal  | 

energy  sot  to  1150  Mb-cc/gm  (Case  l’I5).  The  TOA  \ s depth  curve  I 

is  included  on  figure  ii  and  is  seen  to  fall  f a i r 1 v close  to  the 

^ Ua'eline  Case.  On  ’he  Z vs  Ifj^j  plot  in  figure  15,  Case  015, 


KiKiire  1^-  Kflt'il  ol  luLtriiaJ  KncrKy  of  I'lycr  J’laLes  wil'.i 

\on-miiforni  Vf  loo  i t y Prof  IK-  on  rimo  ol  Arrival. 
(All  Cast's,  V=  -40/70  rm/.jsec  Kamo) 


which  has  a normalized  total  energy  of  1.256,  is  seen  to  have 
the  same  TOA  sensitivity  factor,  E,  as  an  a 1 1 - k i net i c - energy  , 
uniform  velocity  plate  having  a normalized  total  energy  of 
.909.  The  ratio  of  equivalent  kinetic  energy- to  - interna  1 energy 
is  .909/1.250  = .72.  Thus,  internal -energy  in  the  flyer  plate 
is  consistently  equivalent  to  about  3/4  as  much  kinetic  energy, 
insofar  as  the  effect  on  the  shock  introduced  into  the  granite 
block  and  tlie  resulting  TOA  vs  depth  are  concerned. 

4 . 2 Granite  Block  Parameters 

4.2.1  Different  Granite  DOS  Models 

The  Baseline  Case  (BL)  used  the  Pyatt  EOS  model.  In 
Cases  E2  through  E8,  tlie  alternative  EOS  models  described  in 
Appendix  B were  used  for  both  tlie  granite  block  and  tlie  flyer 
plate.  Figure  17  shows  the  Hugoniots  calculated  for  eacli  of  the 
eiglit  EOS  models.  Note  that  the  Hugoniots  for  five  of  the 
models  nearly  coincide  (EASE,  Green-1,  Green-2,  Schuster,  and 
1 i 1 lot  son) . The  Hugoniots  for  the  other  three  models  (Pyatt, 
Green-3,  and  Kalitkin)  arc  somewhat  softer.  Tl.is  behavior 
directly  reflects  in  the  TOA  vs  depth  curves  shown  in  Figure  18; 
the  three  softer  models  (Pyatt;  Case  BE,  Grecn-3:  Case  E4 , 
and  Kalitkin:  Case  E5)  produced  consistently  later  arrival  times 
than  the  five  stiffer  cases.  The  Pyatt  EOS  used  in  the  Baseline 
Case  gives  the  latest  arrival  times  of  all  the  models.  Figure  19 
shows  the  effects  of  the  EOS  models  on  the  plots  of  TOA  sensi- 
tivity factor,  Z,  vs  total  energy,  Substitution  of  alter- 

native models  for  the  Pyatt  EOS  model  has  the  following  effects 
on  Z : 


I 


too. 


E6  (LASL) 


9 

8 


eo. 


— I 1 r 

<H).  eo>  ao. 

TIME  CMICROSEC) 


— r 
100. 


ICO. 


Figure  18.  Time  of  Arrival  Curve  for  Each  of  the  Eight 
EOS  Models  of  Granite.  (50  cm/usec  Uniform 
Velocity  Flyer  Plate). 
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Figure  19.  TOA  Sensitivity  Factor  vs  Normalized  Energy  for  Cases  in 

j whlcli  the  Equations  of  State  of  the  Granite  Block  was  Varied. 
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EOS  Model 

Case 

Effect  on  T 

Equivalent  Change  in 
KE  of  flyer  which  wil 
have  same  effect  on  Z 

Pyatt* 

Basel i ne 

Green- 3 

E4 

.012 

+ 3 a 

Ka 1 i tk  in 

E5 

. 021 

+ 4°. 

Green- 1 

E2 

. 077 

+ 18°„ 

EASE** 

Eti 

.080 

+ 19°o 

Schuster 

E7 

. 004 

+ 24^ 

Green- 2 

E3 

. 094 

+ 24°„ 

Til lotson 

E8 

.114 

+ 30b 

* Used  in  2-D  calculations  by  S-'  and  SAl 
**Used  in  2-D  calculation  by  l.ASI, 

Tlie  maximum  effect  of  tlie  l.OS  model  on  the  TOA  Sensitivity 
factor  (for  the  models  considered)  was  equivalent  to  a 50. 
change  in  the  kinetic  energy  of  the  flyer  plate.  Since  subse- 
c)uent  effects  of  debris  slaji  (e.g.  ground  motions)  are  expected 
to  be  nominally  proportional  to  the  culie-root  of  the  coupled 
energy,  an  uncertainty  in  the  DOS  model  which  is  equivalent 
to  a 30°  uncertainty  in  the  kinetic  energy  of  tiie  debris  is 
jirobably  not  important. 

4.2.2  Variations  within  Pyatt  fOS 

To  examine  the  sensitivity  of  the  TOA  to  variations  in  the 
Pyatt  liOS , four  cases  were  run  in  whicli  the  model  was  perturlied 
or  the  method  of  using  it  was  altered. 

As  outlined  in  Appendix  B,  the  Pyatt  model  for  granite  is 
a complex  formulation  incorporating  the  Salia  theory,  a Birch- 
Murnaghan  fit  of  tlie  zero  pressure  isotherm,  llugoniot  data,  and 
an  ajiproximatc  treatment  of  tlie  120-3.30  kb  qua  rt  z - s t i shov  i t c phase 
change.  To  expedite  its  use  in  code  calculations,  the  model  has 
been  reduced  to  a table  look-up  form.  Rather  than  reconstruct 
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the  table  by  varying  parameters  in  the  basic  formulation,  it 
was  felt  that  the  sensitivity  of  results  to  the  model  could  be 
at  least  crudely  determined  by  arbitrarily  changing  the  pressure 
given  by  the  table  for  the  conditions  in  each  computational 
cell  and  cycle.  In  Case  B1 , this  pressure  was  increased  by  201; 
in  Case  B2,  it  was  decreased  by  20%.  The  effects  on  the  TOA  vs 
depth  curves  in  Figure  20  are  seen  to  be  small.  Reducing  the 
pressure  (Case  B2)  had  a somewhat  larger  effect  on  the  TOA 
tlian  a corresponding  increase  (Case  Bl).  Effects  on  the  TOA 
sensitivity  factor,  Z,  arc  seen  in  Figure  21.  Reducing  the 
pressure  by  20°t,  tuid  an  effect  equivalent  to  reducing  the  flyer 
plate  kinetic  energy  by  4.7d. 

In  another  change  to  the  Pyatt  model,  the  release  adiabats 
in  tlic  model  were  replaced  in  Case  B3  by  a simple  y-law  depen- 
dence i . c . , 

P - (Y,|-npE  (2) 

Y|l  in  Eejn.  2 was  cliosen  for  each  calculat  ional  zone,  based 
on  the  llugoniot  pressure  experincccd.  Thus 

(Yii'l)  = I’h/Ph^H  (5) 

As  seen  in  Figure  22,  the  y'Ipw  release  paths  do  not  reflect 
the  detailed  structure  of  paths  calculated  with  the  Pyatt  model 
for  sliocks  in  the  range  below  about  10  Mb.  .Nonetheless,  tlie 
TOA  vs  depth  curves  for  Case  B3  coincides  almost  exactly  with 
the  Baseline  Case  in  Figure  20. 

The  final  modification  to  the  Pyatt  EOS  model  was  to  change 
the  method  of  interpolation  between  two  values  in  the  tables. 
Normally,  the  logs  of  the  two  quantities  are  first  determined, 
then  linear  interpolations  arc  made  between  logs,  then  the 
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Figure  20.  Variations  in  Time  of  .\rrival  Due  to  Modifications 
of  the  Pyatt  FOS  Model  of  Granite. 
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Figure  22.  Hugoniot  and  Release  Paths  for  Standard  and 
Modified  Pyatt  Granite  Equation  of  State. 


antilog  of  the  resulting  number  is  determined.  In  Case  B6, 
this  jirocedure  was  replaced  by  simple  linear  interpolation, 
directly  between  the  numbers  in  the  table.  The  resulting 
TOA  vs  depth  curve  from  Case  B6  is  virtually  indistinguishable 
from  the  Baseline  Case.  This  suggests  either  that  the  intervals 
between  the  values  in  tlie  Pyatt  EOS  table  are  too  fine  (and  that 
computer  memory  could  be  saved  by  reducing  the  table  size),  or 
tluit  the  costly  process  of  using  logs  and  antilogs  could  be 
replaced  by  simple  direct  interpolation  between  values  in  the 
tables . 

4.2.3  Effects  of  EOS  Model  Under  Otlier  Conditions 

To  test  the  effects  of  the  EOS  model  under  other  conditions, 
comparisons  were  made  between  the  Pyatt  and  EASE  models  for 
cases  in  which  the  velocity  and  ther  internal  energy  level  in 
the  flyer  plate  were  varied.  (The  EASE  EOS  was  selected  for 
tliese  comparisons  because  of  its  use  in  the  EASE  2-D  calculations.) 

The  basic  comparison  of  the  Pyatt  and  EASE  EOS  models  has 


been 

described  in 

Sect  ion  4 , 

,2.1  above . 

In  summary: 

Case 

EOS  Model 

Ve loc  i ty 

Internal 

Energy 

E 

E.qu  i valent 

Kinetic  Energy, 

^•eff 

BE 

Pyatt 

50 

-0 

0 

1 . 0 

E() 

EASE 

50 

-0 

.080 

1 . 19 

Thus 

, use  of  the 

EASE  model 

produces  an 

effect 

on  the  calculated 

TOA  which  is  equivalent  to  a 19°i.  increase  in  flyer  plate  kinetic 
energy . 

When  the  comparison  is  made  at  40  cm/pscc,  the  following 
is  obtained: 
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Case 

EOS  Model 

Velocity 

Internal 

Energy 

Z 

Equivalent 
kinetic  Energy, 

*^eff 

P2 

Pyatt 

40 

' 0 

.212 

.640 

L6B 

LASL 

40 

'0 

.101 

.815 

Under  tliesc  conditions,  the  LASL  EOS  produces  an  effect 
equivalent  to  a 27°o  increase  in  flyer  plate  kinetic  energy. 

When  a -40/70  cm/usec  ramp  velocity  is  used,  together  with 
an  internal  energy  level  of  500  Mb-cc/gin  in  the  flyer  plate, 
these  comparative  results  are  obtained; 


Case 

LOS  Mode  1 

Vcloc  ity 

Internal 

Energy 

Z 

E(]uivalent 
kinetic  Energy 

^eff 

PI  6 

Pyatt 

-40/70 

500 

.015 

1 . 032 

L6A 

LASL 

-40/70 

500 

. 092 

1 .225 

Here,  use  of  the  LASL  LOS  is  equivalent  to  increasing  the  flyer 
plate  kinetic  energy  by  19°o. 

The  above  comparisons  indicate  that  the  effects  of  the 
LOS  model  arc  reasonably  independent  of  the  assumed  velocity 
and  internal  energy  conditions. 

4.2.4  Prclicating  of  Granite  Block 

In  two  cases,  B4  and  B5,  it  was  assumed  that  the  granite 
block  is  preheated  by  direct  radiation  from  the  nuclear  source. 
F-or  these  runs,  the  internal  energy  density,  lE^,  at  the  front 
surface  of  the  block  was  specified  to  be  0.2  and  20  Mb-cc/gm. 
Based  on  results  of  a preliminary  S'^  analyses,  it  was  further 
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specified  tliat  the  internal  energy  density  in  the  block  falls 
off  witn  depth,  D,  within  the  block  according  to 


Ilin  = IE  X 
D o 

The  total  initial  internal  energy  added  in  the  granite  block 
due  to  preheat  in  Case  B4  was  approximately  equal  to  the 
internal  energy  added  to  the  flyer  plate  in  Case  Pll.  Similarly, 
the  added  internal  energy  in  the  block  in  Case  B5  was  about  the 
same  as  that  added  to  the  plate  in  Case  P12.  (Cases  Pll  and  P12 
were  discussed  in  Section  4.1.3.)  It  is  therefore  very  inter- 
esting that  the  TOA  vs  depth  curves  for  Cases  B4  and  Pll,  and  for 
Cases  B5  and  P12,  almost  exactly  coincide,  as  seen  in  Figure  14, 
on  page  29.  The  Z vs  E^^^  values  for  these  cases  also  coincide 
(Figure  15,  on  page  30).  Thus  adding  internal  energy  to  the 
granite  block  has  the  same  effect  as  adding  the  same  internal 
energy  to  the  flyer  plates.  In  either  case,  adding  a given 
amount  of  internal  energy  is  the  equivalent  of  adding  about  3/4 
of  the  same  amount  of  kinetic  energy.  This  again  indicates 
that  while  the  total  amount  of  energy  available  has  an  impor- 
tant effect  on  the  TOA  vs  depth  in  the  granite  block,  the 
specific  distribution  of  tlie  energy  has  only  a minor  effect. 

4 . 3 Numerical  Parameters 

4.3.1  Zone  S i zc 

All  of  the  preceding  cases  used  1-D  cells  with  initial 
thicknesses  of  0.4  cm  in  the  flyer  plate  and  0.4  cm  near  the  face 
of  the  granite  block.  The  tliickness  of  cells  in  the  block 
increased  in  0.5°^  increments.  These  thicknesses  were  approxi- 
mately the  same  as  had  been  used  in  the  2-D  analyses  by  S^, 

LASL,  and  SAI.  To  test  the  adequacy  of  this  zoning.  Case  Cl 
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was  run  in  which  all  initial  thicknesses  were  cut  in  half. 

The  effect  on  TOA  vs  depth  was  quite  small:  cutting  the  zone 
size  in  half  had  an  effect  equivalent  to  increasing  the  kinetic 
energy  by  2.21.  This  suggest  that  the  zoning  in  the  2-D  studies 
was  adequate. 

4.3.2  Artificial  Viscosity 

Artificial  viscosity  is  used  in  finite  difference  code 
solutions  to  provide  stability  across  the  discontinuity  of  the 
shock  front.  It  has  the  effect  of  smearing  the  front  over 
a few  cells,  and  tends  to  reduce  the  peak  pressure.  In  Case  C2, 
the  artificial  viscosity  parameter,  Q,  was  changed  from  4 (the 
value  used  throughout  this  study)  to  1.  The  effect  on  TOA  vs 
depth  was  negligible  - equivalent  to  increasing  the  kinetic 
energy  of  the  flyer  plate  by  about  1% . 

4.3.3  Energy  Iteration 

Since  all  of  the  codes  used  in  the  various  HUSKY  PUP 
studies  were  explicit  in  tlieir  formulation,  there  is  a problem 
in  calculating  both  pressure  and  energy  at  the  same  time.  Many 
codes  approach  this  problem  by  iterating  the  calculations  of  P 
and  c until  they  converge.  In  CRALE-1,  the  final  pressure  and 
its  derivative  with  respect  to  e are  approximated  by  calculating 
a value  based  on  the  energy  from  the  previous  cycle  and  then 
adding  the  incremental  change  produced  by  the  change  in  energy 
(-Pd\').  When  this  correction  term  was  removed  in  Case  C3,  the 
effect  on  the  TOA  vs  deptli  was  negligible,  suggesting  that  the 
potential  error  which  could  result  from  improper  time-centering 
of  the  energy  used  in  the  pressure  calculation  is  quite  small. 
The  iterative  procedures  used  in  some  codes  can  be  expensive, 
particularly  if  the  EO.S  model  is  very  complex;  the  results  of 
Case  C3  suggest  these  procedures  could  be  simplified  or  even 
eliminated  without  serious  loss  of  solution  accuracy. 

J8 


Figure  23  compares  all  of  the  TOA  vs  depth  curves  obtained 
when  the  above  changes  in  the  numerical  parameters  were  made 
(Cases  Cl,  C2,  and  C3) . The  effects  are  almost  indescernible , 
indicating  that  these  aspects  of  the  numerics  are  not  likely 
to  have  significantly  influenced  the  2-D  calculations  of  TOA 
in  the  granite  block. 
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Figure  2i.  Variations  in  Time  of  Arrival  Due  to  Cl-anges 
in  Numerical  Techniques. 


SECTION  5 


DIFFERENCES  IN  EQUIVALENT  KINETIC  ENERGY  BETWEEN 
EXPERIMENTAL  DATA  AND  2-D  CALCULATIONS 

It  is  of  interest  to  determine  the  TOA  sensitivity  factor, 
Z,  and  equivalent  effective  kinetic  energy,  E^j-j-,  for  the 
experimental  TOA  data  and  for  the  values  obtained  from  the  2-D 
calculations.  This  is  done  by  first  calculating  E by  comparing 
these  TOA  vs  depth  curves  with  the  curve  from  the  Baseline  Case 
(Figure  3 on  page  10  and  Figure  6 on  page  15).  Inserting  these 
values  for  E in  Figure  24,  it  is  then  possible  to  estimate  the 
kinetic  energy  of  the  1-D  flyer  plate  which  would  produce 
essentially  the  same  TOA  vs  depth  curves  as  are  seen  in  the 
experiment  and  in  the  results  of  the  2-D  calculations.  The 
absolute  value  of  this  equivalent  kinetic  energy  is  not  of 
itself  very  significant,  but  the  relative  differences  are  a 
measure  of  how  close  the  2-D  and  1-D  calculations  came  to  the 
measured  data. 


TOA 

Sens i t i v i ty 
Factor , 

E 


Basel  ine  1 -1)  Case 

1.0 

HUSKY  PUP 

Measured  Data 

- .047 

LASL  2-D 

Calculation  (pre-test) 

- . 135 

S^  2-D 

Calculatidn  (pre-test) 

- . 174 

SAI  2-D 

Calculat ion (pre - test ) 

- . 46 

Equivalent  Difference  of 

Effective  KE,  Eeff  from  Value 
Ep£-j-  for  Measured  Data 

1.0  10?i 

.91 

.75  181 

.70  21% 

.33  641 


.51 


i 

The  LASL  and  S-'*  2-1)  results  are  probably  within  the 
uncertainty  limits  of  the  HUSKY  PUP  source  yield.  Tiiis  fact, 
plus  the  findings  of  the  current  study  that  calculations  of 
■fOA  vs  depth  curves  are  insentitive  to  relatively  large  changes 
in  input  models  and  assumptions,  indicates  that  further  refine- 
meTits  of  the  LASL  and  S^  2-0  calculations  arc  not  necessary 
nor  likely  to  be  fruitful.  We  conclude  that  the  tcclmology 
demonstrated  by  these  2-D  calculations  is  addciuatc  to  jircdict 
debris  slap  for  cratering  and  ground  motion  analyses. 

The  major  difference  between  the  SAl  2-0  results  and  the 
measurements  and  the  results  of  the  otlier  2-0  calculations 
are  not  exp  1 a i nali  1 c , insofar  as  we  can  tell,  by  differences  in 
the  input  assumptions.  Notic  of  the  variables  or  uncertainties 
wliich  we  considered  in  the  current  study  appear  likely  to  be 
responsible  for  the  large  discrepancy  in  the  SAl  results. 

This  suggests  the  i^ossibility  of  a ca  Iculat  iona  1 error  of 
some  type. 


55 


APPENDIX  A 


A SIMPLE  INTERPRETATION  OP  THE  HUSKY  PUP  TIME-OF-ARRIVAL  DATA 

Attempts  were  made  in  HUSKY  PUP  to  measure  pressures  and 
velocities  at  various  depths  in  the  granite  block.  However, 
the  only  dependable  data  obtained  were  the  t imes-of-arrival 
(TOA)  of  the  first  shock  at  the  various  gages.  These  data, 
shown  in  Figure  Al,  are  quite  consistent;  a least  squares  fit 
of  all  the  data  in  the  form 


D = at^''  (Al) 

has  an  RMS  error  of  less  than  2°6  , independent  of  whether  depth, 
1),  is  measured  from  the  front  surface  of  the  block  or  the  center 
of  tlie  source.  The  error  drops  to  about  1%  if  only  the  Physics 
International  (PI)  ladder  data  are  fit,  and  is  smaller  vet 
if  either  the  earliest  points  arc  excluded,  or  an  additive 
constant  is  used,  i.e.,  1)  = at^’  + c.  The  coefficients  a and 
b fall  into  two  classes,  depending  on  whether  the  distance  to 
the  datum  point  is  measured  from  the  front  surface  of  the  block 
or  the  center  of  the  device.  Table  Al . 

•A  pressure  vs  range  relationship  can  be  derived  by  com- 
bining the  fits  to  the  TOA  data  with  Hugoniot  shock  velocities, 
U^,  obtained  from  the  EOS  models.  Curves  of  U^,  vs  P are  shown 
in  Figure  A2  for  two  of  tlic  models.  The  shock  velocities  at 
anv'  pressure  above  1 Mb  differ  by  less  than  5%  between  the  eight 
EOS  models  used  in  this  study  and  can  be  approximated  in  the 
range  .7)  to  100  Mb  by  a power  law,  namely 


Us  = dP*^ 


(A2) 


where  d = .89  and  e = .44.  This  fit  is  included  in  Figure  A2. 
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TABLE  A1 


ALTERNATIVE  LEAST  SQUARES  FITS  TO  HUSKY  PUP 
TOA  DATA  AND  INFERRED  PRESSURE-DEPTH  RELATIONSHIPS 

(D  = at*’  + c)  (TOA  data) 

Q 

(P  = aCD-c)*^  (inferred  P-D  relationship) 


I.  Where  D is  measured  from  center  of  device: 


a 

b 

c 

a 8 

Using  all  27  data  points; 

36.18^^^ 

.383^*^^ 

o(i) 

8*  * 
2.61x10®  -3.66 

or ; 

16.48 

.517 

35.31 

6.5  xlO^  -2.12 

Using  only  PI  ladder  data: 

36.99 

.376 

0 

or : 

13.95 

.553 

40.03 

Using  only  last  7 ladder 

points : 

34.86 

.392 

0 

or: 

25.25 

. 446 

16.97 

II.  Where  D is  measured  from 

front  surface  of 

granite 

block : 

Using  all  27  data  points: 

19.20^^^ 

.490^^^ 

0(2) 

** 

2.3  X 10^  -2.366** 

Usingly  only  PI  ladder  data: 

19.00 

.494 

0 

or : 

15.25 

.536 

7.34 

Using  only  last  7 ladder 

points : 

18.67 

.499 

0 

or : 

18.66 

.499 

.004 

(1)  curve  "Y"  on  Figure  A1 

(2)  curve  "Z"  on  Figure  A1 

* curve  P]^  on  Figure  A3 

**  curve  P2  on  Figure  A3 
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'd-cI 

a 


(A3) 


whore  a,  li , c are  the  coefficients  listed  in  Table  Al. 

can  be  eliminated  in  liquations  A2  and  A5  above  to  produce  a 

pressure- d i stance  rel  at  ionsh  i of  the  form 


P = udl-c)*^  (A4) 

where  a = (b/d)a^'^*’  and  3 = (!^ — i-)  (— ) 

b e 


Two  of  tiie  pressure  vs  deptli  curves  derived  from  least  squares 
fits  of  all  tlie  data  points  are  presented  in  Figure  A3.  Since 
tl\e  curves  represent  two  different  fits  to  the  same  data,  it 
IS  not  surprising  tliat  they  cross  in  the  middle  of  tlie  region 
of  data.  llie  small  difl'erence  between  shock  x'elocities  in  tlie 
\ai-ious  laTS  models,  jilus  the  rather  small  difference  between 
the  pressure  vs  depth  fits  in  ligui'e  A3  (less  than  30).  for 
all  pressures  above  (1.7  Mb),  suggest  tliat  the  higli  pressure 
I’-D  relation  is  relati\ely  insensitive  to  the  HOS  model  or  to 
the  accuracy  of  individual  data  points. 


APPENDIX  B* 


Tin;  GRANITE  EQLIATION-OE-STATE  FOR  TEMPERATURES 
BETWEEN  0.5  EV  and  1 K'EV 

B. 1 Introduct ion 

Eight  models  of  the  quartz  (Si02)  equa t i on  - of - stat e (EOS) 
have  been  studied  to  estimate  tlic  potential  error  in  prediction 
calculations  of  tlie  IIUSKV  PUP  granite  block  experiments  caused 
by  uncertainties  in  high  energy  material  properties.  Of 
specific  concern  was  tlie  material  behavior  in  tlie  temperature 
range  between  0.5  e\'  and  1.0  keV.  Radiation  effects  dominate 
the  calculation  above  1 keV  (10^-'’  ergs/gm),  while  the  various 
experiments  of  interest  in  IIUSKV  PUP  occur  at  temperatures  above 
0.5  e\'.  Tbe  temperature  range  of  interest  is  well  above  that 
of  most  previous  experiments,  hut  it  is  too  low  to  expect  simple 
Thomas- Perm i or  Saba  models  to  be  entirely  valid.  Of  tbe  eight 
specific  models  investigated  (fable  B1  and  References  5 through 
12), one  was  based  on  tbe  Saba  equation  (5);  tbe  other  seven 
used  tbe  Thomas  Fermi  theory  (12)  with  various  corrections  and/or 
approximations.  Comparisons  of  the  Hugoniots  and  release  ad  i - 
abats  for  all  of  tbe  models  arc  presented  in  Section  B.  2 , while 
specific  clia  rac  ter  i St  ics  of  each  model  arc  discussed  briefly 
in  Sect  ion  B . 5 . 


* This  appendix  was  originally  a memo  written  while  one  of  tbe 
authors  (Schuster)  was  at  R 5 D Associates  and  represents 
work  sponsored  there  by  DN'A . 
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Tabic  Bl.  Si02*  iiOS  MODliLS  STDDII-.D 


AUTHOR 

1 

THEORY 

TYPE 

REPERENCE 

1.  PYATT  (S^) 

SAHA 

TABLE  LOOKUP 

5 

2.  BARYNH  (LAST) 

THOMAS- PERM I 

1 1 

0 

5.  GRPEN  1 (RDA) 

ft 

11 

7 

4.  " 2 (RDA) 

" a 

1 1 

- 

5.  '•  3 (RDA) 

" b 

1 1 

•7 

6.  KALITKIN 
(USSR/LLl) 

" c 

1 1 

8.9 

7.  SCHUSTER  (RDA) 

1 1 

ANALYTIC  PIT 

10 

8.  TIELOTSON  (GA) 

1 1 

.. 

1 I 

*Pq  = 2.b5  gni/cin' 


‘Includes  nuclear  correction  assuming  nuclei  behave  as  a perfect 
gas. 

^Includes  nuclear  correction  assuming  nuclei  behave  as  harmonic 
osc  i 1 lators . 

^■Includes  electron  exchange  and  correlation  corrections. 


B.2  Comparisons  of  Various  Models 

Since  dissipation  from  eneigy  densities  lietween  It)**  and 
10^^  ergs/gm  always  involve  shock  processes,  the  principal 
llugoniot  is  a useful  curve  for  comparing  various  models, 
ihe  jum)!  conditions  at  a shock  f ront , comb  ined  with  eacli  KOS 
model,  arc  sufficient  to  determine  the  locus  of  points  attainable 
in  a strong  shock,  i.e.,  the  principal  llugoniot.  The  jump  con- 
ditions arc: 


P(S-II) 

(conservat  ion 

0 f 

mass) 

(Bl) 

U“/2 

(consc rvat  ion 

0 f 

energy) 

(B2) 

PoSl) 

(.conservat  ion 

of 

momentum) 

(R3) 

where  and  P are  the  pre-  and  post-shock  densities  of  the 
material,  U and  S arc  the  material  and  shock  speeds,  E is  the 
energy  density,  and  P is  the  shock  pressure.  Strictly  speaking, 
P,  E,  and  11  rci^rcscnt  changes  in  those  quantities,  but  along 
the  principal  llugoniot  the  initial  jircssure,  energy  and  mater- 
ial velocity  of  the  rock  arc  small  and  can  be  ignored.  .Since 
the  above  three  ccpiations  contain  fi  vc unknowns  , adding  a fourth 
equation,  namely  the  EOS  model,  allows  one  to  solve  for  any 
four  of  the  unknowns  as  a function  of  the  fifth.  fhe  resulting 
llugoniots,  represented  as  P vs  p,  arc  shown  in  figure  Hi. 

All  of  the  curves  coalesce  near  1 mbar  since  they  were  forced 
to  match  tiic  low  pressure  (<.S  Mb)  experimental  data. 

At  higher  pressures  tiic  Kalitkin,  Pyatt,  and  Orcen  S models  arc 
considerably  softer,  i.e.,  more  compressible,  than  the  other 
five  models  which  remain  in  ;i  tight  band.  At  very  high  pressures 
all  of  the  models  approach  an  effective  gamma  [y^  = (P  pE)/pE] 
of  -S/5  except  for  the  Tillotson  model  which  converges  to  a 
value  of  3/2.  While  there  appear  to  be  rather  striking  differ- 
ences between  the  various  llugoniots,  the  energy  deposited  in  the 


J 


(Mbar) 


material  at  any  given  pressure  does  not  vary  greatly  even  be- 
tween the  extremes  of  the  curves.  Expressing  the  llugoniot 
energy  as  a fun>.tion  of  pressure  and  effective  gamma,  namely 


E 


P 

P (y  + 1 ) 
o ^ ' e 


(B4) 


it  follows  that  at  any  specific  value  of  the  pressure,  the 
ratio  of  deposited  energies  between  any  two  models  is  simply 


R(E) 


Y2M 


(B5) 


Although  in  tneory,  can  vary  between  1 and  for  pressures 
above  5 mbars  tlie  effective  gamm;  along  the  llugoniot  for  any 
of  the  models  studied  was  never  less  than  1.44  or  greater  than 
2.6*.  In  addition,  at  any  pressure,  the  effective  gammas  for 
the  extreme  llugoniots  never  differ  by  more  than  4()o,  so  that 
the  maximum  variation  in  energy  is  less  than  12o,  Figure  B2. 

The  ratio  of  shock  velocities  is  even  less  sensitive  to  differ- 
ences in  the  llugoniots,  varying  as  VR(E).  I'hus , the  differences 
between  the  mod.'is  shown  in  Figure  B1  are  probably  not  impor- 
tant to  a calculation  of  the  overall  material  response,  parti- 
cularly if  one  considers  that  the  llugoniot  is  only  the  loci  ot 
states,  only  one  of  which  is  reached  by  any  specific  zone  of 
material  in  a calculation.  Of  much  greater  potential  signifi- 
cance are  any  systematic  differences  in  release  adiabats  from 
the  various  llugoniots.  While  the  energy  deposited  by  a given 
magnitude  shock  does  not  vary  greatly  between  the  models,  if 


Below  5 mhar  the  material  is  a solid  and  cannot  be  realistically 
represented  by  a y-law  equation. 
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the  fraction  of  the  energy  recovered  during  the  adiabatic  expan- 
sio'i  was  to  differ,  the  effect  could  be  dramatic.  Such  does  not 
appear  to  he  the  case,  as  release  adialiats  from  10,  100,  and 
1000  Mbars  (l-'igure  B3)  , appear  to  cluster  tightly  in  I’-p  space, 
f^ven  when  plotted  on  the  energy- vo  lume  plane  (Figure  B4),  the 
differences  for  release  from  100  and  1000  Mbar  arc  not  great. 

The  energies  lost  in  waste  heat  after  expansion  to  100  cm-^/gm 
differ  by  a factor  of  4 between  the  extremes  of  the  models. 
Ilowevei',  tliese  energies  arc  an  order  of  magnitude  less  than  tiic 
initial  energy  deposited,  so  the  enerjy  reaouered  only  varies 
between  01  and  07.5".  I-'u  r t he  rmo  re  , the  spread  between  the 
presumably  more  accurate  six  table  look-up  models  is  only  2°o 
(0  5 . 5 - 0 7 . 5 0 ) . Tlie  larger  spread  for  tlie  10  Mbar  unloading  is 
due  to  tlie  differences  between  tlie  models'  treatment  of  the  solid 
pliase  of  SiO-)  and  while  important  for  low  stress  attenuation,  was 
not  a subject  of  this  review. 

An  additional  set  of  unloading  paths  from  the  jioint  , = 

.^^tl.t'5  g/cm-^)  and  F,  = lu'^  ergs/ cm,  shown  in  Figure  Bi,  simu- 
lates release  after  the  instantaneous  deposition  of  energ\’  into 
cold  material.  While  not  an  exact  reiiresentat  i on  of  radiation 
deposition,  this  state  is  tyjiical  of  the  initial  conditions 
fretiuentl)'  used  in  calculations  ol  nuclear  explosions.  I'he  six 
table  look-up  moilels  return  between  88  and  92  percent  of  the 
initial  energ}-  during  expansion  to  llld  cm-'’/gm,  while  the  two 
analytic  fits  1 1' i 1 1 ot  son  ' s and  Schuster's),  recover  '.14  and  97. 5'’,., 
I'espec  t ive  1 >■ . If  one  assumes  that  the  energy  recovered  and 
available  t o do  cont i nu i ng  work  is  the  most  significant  driving 
mechanism  of  the  ground  motions,  the  differences  exhibited  in 
Figures  B1  through  B4  do  not  apjiear  significantly  large. 

Two  further  tiucstions  remain:  First,  do  the  models  pre- 
sented here  reall\-  cover  the  complete  range  of  possible  mater- 
ial behavior,  i.e.,  could  the  real  material  load  and/or  unload 
along  juiths  outside  the  limits  of  our  models?  Second,  how 
critical  is  the  time  scale  over  which  the  material  unloads? 

o' 
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P (Mbar) 


E(ergs/gni) 
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Figure  BA.  Release  Adiabats  from  -v-lO,  100  and  lOOOi'B’ar  Hugoniot 

1 5 

Pressures  and  from  p = p^,  E = 10  ergs/gm  for  Quartz 
(Several  curves  arc  deleted  where  they  coalease  with 
the  others) 

(>;) 


Clearly,  if  tlic  unloading  wave  speeds  of  one  model  are  signi- 
ficantly faster  tlian  another,  the  energy  would  lie  released 
quicker  and  could  change  the  motion  at  lower  jieak  jiressures. 

I'he  question  of  whether  one  has  covered  tlie  complete  set 
of  possible  material  models  is  difficult  if  not  impossible  to 
answer;  liowever,  the  issue  lias  been  addressed.  One  simple 
way  to  look  at  the  models  is  to  express  the  et|uat  i on  - of  - sta  te 
in  term^  . f an  effective  gamma,  i'hus , in  figure  H5,  vs 

p for  unloading  from  p = p^ , 1;  = 10*^  and  from  the  llugoniot 
pressure  of  about  1000  Mbar,  is  plotted  for  the  two  extremes 
of  the  table  lookup  models  (Green's  tables  of  the  Thomas- I’ermi 
model  with  the  perfect  gas  and  harmonic  oscillator  nuclear 
corrections)  and  the  two  analytical  fits.  y^  varies  between 
1.2  and  1.7  over  the  pressure  range  0.02  to  1000  Mbar  for  the 
release  jiatlis  iilotted.  While  gamma  can  theoret  tea  1 1 >'  vary 
from  1 to  O’,  it  is  difficult  to  imagine  a material  whose  gamma 
is  consistently  below  1 . .i  or  above  1.7  for  the  energy  densities 
of  interest.  'fhus,  tlie  static  analysis  of  tiie  SiO,  equation- 
of-state  models  presented  above  strongly  suggests  that  for 
temiieratures  above  0..S  e\',  details  of  tlie  specific  COS  used 
piohably  will  not  be  significant  in  the  calculations  of  material 
response  to  a nuclear  explosion. 

H.  .1  Individual  Models 

I'he  eight  models  for  quartz  compared  in  the  previous  section 
are  discussed  individually  below.  While  a brief  description  of 
the  aspects  of  each  is  presented,  this  is  not  intended  to  he  an 
in-depth  review.  .All  of  the  models  assume  density  and  specific 
energy  are  the  indeiiendent  variables.  More  comprehensive  dis- 
cussions of  each  model  may  be  found  in  the  appropriate  references 
listed  in  Tab  1 e lU  . 


B. 3. 1 Pyatt  Model  (S^) 


The  CAPO  code,  used  to  generate  the  basic  data  from  which 
these  tables  were  derived,  solves  the  Saha  equation  assuming 
local  thermal  equilibirum  and  includes  electron  excliange  and 
correlation  terms.  The  low  temperature  behavior  is  presented 
by  a Birch-Murnaghan  type  equation  with  parameters  chosen  to 
match  experimental  Hugoniot  data  up  to  120  kbar.  The  quartz- 
stishovite  phase  change  which  occurs  at  shock  pressures  between 
120  kbar  and  350  kbar  was  included  at  the  120  kbar  level,  so 
the  calculated  principal  Hugoniot  is  consistently  lower  than  the 
measured  data  in  tliis  region. 

I’ressure  and  temperature  are  each  represented  by  a coarse 
and  fine  table  in  tliis  model.  This  two  table  scheme,  in  which 
the  fine  table  is  embedded  in  the  coarse, was  necessary  to 
represent  the  phase  change  without  a prohibitively  large  com- 
puter storage  requirement.  The  coarse  tables  were  dimensioned 
23  X 38  (density  x energy)  and  included  densities  between  10'^ 

! and  25  gm/cm-^  and  energies  between  1.58  x 10^  and  4 x 10^^  ergs/gm. 

j The  fine  table  (53  x 53)  included  densities  between  1 and 

I 10  gm/cm^  and  energies  from  10^  and  10^^  ergs/gm.  Thus  only 

17400  entries  were  needed  for  both  pressure  and  temperature. 

(A  comparable  single  table  model  generated  earlier  needed  23,000 
entires  to  provide  tlie  same  coverage  and  was  too  largo  for 
several  of  the  codes  to  use  conveniently.)  A linear  interpolation 
using  the  logarithms  of  all  the  entries  was  employed  to  obtain 
pressure  and  temperature  for  densities  and  energies  within 

(the  bounds  of  the  tables.  Outside  the  table  limits  the  pressures 
were  calculated  by  extrapolating  from  tlie  nearest  table  entries 
assuming  a perfect  gas  behavior.  (This  same  procedure  was  used 
in  all  the  table  look-up  models  presented  here.) 


B.3.2  Barnes  Model  (LASL) 
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The  standard  LASL  EOS  is  derived  from  the  Thomas-Fermi  (T-F) 
formulation  for  an  electron  gas.  The  solid  behavior,  high- 
temperature  nuclear  motion  and  electron  exchange  and  correla- 
tions forces  are  included.  F’ressure  and  temperature  are  found 
by  a two-step  process  in  this  model.  First,  the  pressure,  P^, 
and  energy,  E^ , at  the  desired  density,  are  calculated  for  the 
zero  degree  isotherm  by  interpolating  in  two  33  entry  tables. 

The  thermal  contribution  to  pressure  and  the  temperature  are 
then  calculated  from  33  x 16  tables  (density  x energy)  using 
the  density  and  thermal  energy  (E-E^).  The  tables  include 
densities  and  energies  from  0.02  to  14.4  gm/cm"^  and  0 to  2.64 
X 10^“^  ergs/gm  respectively.  Since  the  tables  for  thermal 
energy  start  from  zero,  a linear  interpolation  of  the  actual 
values  rather  than  their  logs  was  necessary. 

B.3.3  Green  1 Model  (RDA) 

• J.  Green  (RDA)  derived  data  for  tiiis  model  by  fixing 
the  temperature  and  the  electron- free  energy,  and  using  Latter's 
T-F  results  to  obtain  the  resulting  pressure,  energy,  partial 
volumes,  and  entropy  contributed  by  the  electrons  from  both  the 
silicon  and  the  oxygen  atoms.  Since  the  electrons  of  both  atoms 
are  in  thermodynamic  equilibrium  if  they  share  a common  temper- 
ature and  electron- free  energy,  the  electronic  pressures  are 
the  same  at  the  atomic  interfaces.  The  pressure  of  tlie  mixture, 
therefore,  is  the  pressure  at  the  edge  of  either  atom. 

The  original  density,  pressure, and  energy  output  was  con- 
verted into  30  X 27  tables  (p  x E)  of  pressure  and  temperature 
by  an  auxilliary  routine.  In  addition,  since  the  T-F  model 
used  is  only  valid  for  a hot  electron  gas,  at  energy  densities 
less  than  3 x 10^1  ergs/gm  the  pressure  was  calculated  using  the 
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low  pressure  portion  of  the  Schuster  analytic  fit.  A simple 
weighting  function  was  used  to  go  smoothly  between  the  fit 
and  tlie  tallies  for  energies  lietween  0.3  and  5.0  x 10  ergs/gm, 

i . e.  . 

■ f’.-it  * “•f’table  O'-) 


where 


f 


- 0-3  X 10^' 
3 X 10l2 


(B7) 


I'he  logarithmic  interpolation  routine  used  for  the  I'yatt  tables 
was  also  used  for  all  the  (ireen  ^ Kalitkin  tables. 


B.5.4  Green  2 Model  (kllA) 

Ihe  second  Green  model  differs  from  the  one  .■ibo\e  in 
that  the  effects  of  nuclear  motion  are  included  by  assuming 
the  nuclei  behave  as  a jierfect  gas  with  three  degrees  of  freedom. 
Ibis  results  in  an  increase  in  both  pressure  and  energy  at  a 
given  temperatuie  and  shifts  the  llugoniot  slightly  to  the  left 
(less  compressible)  in  I'-p  space  (figure  Bl). 

B.3..3  Green- 3 Model  (RDA) 

Again,  the  original  Green  model  was  motlified  by  a nuclear 
motion  contribution,  this  time  assuming  the  nuclei  are  bound  by 
a harmonic  well.  The  nuclei,  therefore,  contribute  nothing 
to  the  pressure  but  add  3 kt  per  nucleus  to  the  energy,  pro- 
ducing a significant  softening  of  the  llugoniot,  shifting  it 
to  the  right  in  figure  Bl. 

The  two  models  of  nuclcai-  motions  aie  exceedingly  useful 
since  they  should  represent  ui^per  and  lower  limits  of  this 
effect  on  the  pre>sure.  As  shown  in  figure  Bl,  above  .30  Mbar 
the  llugoniots  produced  by  the  Cirecn- 2 and  -3  models  do  indeed 
bound  the  other  curves. 
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B.3.b  Kalitkin  Model 
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Corrections  to  the  basic  T-I-  description  for  tlie  electron 
quantum  and  exchange  forces  at  ()°K  were  derived  hy  Kalitkin 
These  clianges  were  extended  at  LLL  and  incorporated  in  a code 
which  generates  tables  of  pressure,  energy,  and  correction 
terms,  AP  and  AF.  as  functions  of  densit)’,  temperature  and  atomic 
number.  By  assuming  an  average  atomic  weight  of  10.662,  appro- 
priate liigh  energy  tables  were  generated  for  quartz.  Again  the 
Schuster  fit  was  used  at  low  enei'gies. 

The  average  atom  ajiprox  imat  ion  is  cruder  tlian  equating 
temperature  and  free  energy  as  in  the  Creen  models.  However, 
setting  AP  and  AF  to  zero  i)roduced  an  FOS  which  overlays  that 
of  Creen-l,  so  t lie  two  methods  ajqiear  et|uivalent  for  quartz. 

Tlie  Kalitkin  model,  therefore,  can  he  comjiared  with  the  other  T-F 
models  to  estimate  the  effects  of  the  electron  correction  terms. 

B..3.“  Schuster  Model  (RDA) 

In  an  attempt  to  get  a simple  anal\-tic  FOS  for  ciuartz,  an 
expression  was  found  to  a]!]!  I'ox  i ma  t e the  effective  gamma  calcu- 
lated from  the  Cireen-l  model.  The  resulting  exiiression 

(r-l)  = (0..3.S  log  - - 0.4641“  + 0.4  + 0.12  log  p (B8) 

P 

jirovides  a gas  iiressure  lassuming  a iierfect  gas)  which  is  then 
.idded  to  the  solid  phase  contribution.  The  resulting  llugoniot 
and  associated  release  adiahats,  Figures  B1  and  B.3,  are  in 
reasonable  agreement  with  those  of  Oreen-l.  The  unloading  1 rom 
normal  density  and  1.  = lo'-''  ergs/gm  differ  markedly,  however, 
suggesting  a better  fit  should  be  found  if  this  model  is  to  be 
used  in  future  calculations. 
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B.3.8  Tillotson  (GA)  Model 


An  analytical  fit,  first  derived  by  J.  Tillotson  for 
hypervelocity  impact  studies  of  metals  has  been  used  extensively 
in  calculations  of  motions  in  rocks.  The  fit  has  two  parts,  one 
which  models  the  solid  behavior  (using  a polynominal)  and  the 
second  which  represents  the  gasous  state.  The  gas  term  has  a 
perfect  gas  form  P = (y-13pE>  in  which  the  gamma  is  allowed  to 
vary  thusly: 


(Y-1) 


a + 


EoH' 


r+  1 


(B9) 


At  high  energies  this  equation  reduces  to  yl  = . where  a has 

historically  been  assumed  to  be  0.5.  As  seen  in  (•igurc  Bl, 
the  resulting  Hugoniot  is  consistent  with  several  of  the  other 
models  up  to  1000  Mbar.  Above  this  pressure  the  crfect.>''e 
gamma  of  the  other  models  begin  to  increase  asymptotically  to 
5/3  while  the  Tillotson  value  remains  at  1.5 
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ATTN;  J.  G.  Lewis 
ATTN:  Harold  L.  Brode 
ATTN-  Cyrus  P.  Knowles 
ATTN:  William  B.  Wright,  Jr. 

ATTN:  Henry  Cooper 
ATTN:  Jerry  Carpenter 
ATTN;  Albert  L.  Latter 

Sciences  Applications,  Inc. 

ATTN:  David  Bernstein 
ATTN:  D.  E.  Maxwell 

Science  Applications,  Inc. 

ATTN:  Tech.  Library 

Southwest  Research  Institute 
ATTN:  Wildred  E.  Baker 
ATTN:  A.  B.  Wenzel 

SRI  International 

ATTN:  Burt  R.  Gasten 
ATTN:  George  R.  Abrahamson 
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Systems,  Science  and  Software,  tnc. 

Universal  Analytics,  Inc. 

ATTN 

Tech.  Lib. 

ATTN:  E.  I.  Field 

ATTN 

Donald  R.  Grine 

ATTN 

Ted  Cherry 

Weidlinger  Assoc.  Consulting 

Engineers 

ATTN 

Thomas  D.  Riney 

ATTN:  J.  W.  Wright 

ATTN:  Melvin  L.  Baron 

Terra  Tek, 

Inc. 

ATTN 

Tech.  Lib. 

Weidl inger  Assoc.  Consulting 

Engineers 

ATTN 

Sidney  Green 

ATTN:  J.  Isenberg 

Tetra  Tech 

, Inc. 

Westinghouse  Electric  Corp. 

ATTN 

Tech.  Lib. 

Marine  Division 

ATTN 

Li -San  Hwang 

ATTN:  W.  A.  Volz 

TRW  Defense  & Space  System  Group 

ATTN 

D.  H.  Baer,  Rl-2136 

2 cy  ATTN 

Peter  K.  Dai,  Rl-2170 

ATTN 

I.  E.  Alber,  Rl-1008 

ATTN 

Tech.  Info.  Center,  S-1930 

ATTN 

R.  K.  Plebuch.  Rl-2078 

TRW  Defense  A Space  Sys.  Group 
San  Bernardino  Operations 

ATTN:  E.  Y.  Wong,  527/712 


